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Dimerized quantum spin systems may appear under several circumstances, e.g by a modulation 
of the antiferromagnetic exchange coupling in space, or in frustrated quantum antiferromagnets. In 
general, such systems display a quantum phase transition to a Neel state as a function of a suitable 
coupling constant. We present here two path- integral formulations appropriate for spin 5 = 1/2 
dimerized systems. The first one deals with a description of the dimers degrees of freedom in an 
S0(4) manifold, while the second one provides a path-integral for the bond-operators introduced 
by Sachdev and Bhatt. The path-integral quantization is performed using the Faddeev-Jackiw 
symplectic formalism for constrained systems, such that the measures and constraints that result 
from the algebra of the operators is provided in both cases. As an example we consider a spin- 
Peierls chain, and show how to arrive at the corresponding field-theory, starting with both an S0(4) 
formulation and bond-operators. 
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I. INTRODUCTION 



Dimerized quantum spin systems appear frequently in quantum antiferromagnets, leading to a 
spin-liquid in the form of a valence-bond solid (VBS) for strong enough dimerization. A well known 
example is given by spin-Peierls systems, where in quasi one-dimensional materials the dimerization 
can arise due to electron-phonon interaction [l|, [2f . Dimerized systems constitute also a starting 
point to study frustrated quantum antiferromagnets in two dimensions. In such a way, models 
with antiferromagnetic exchange interactions beyond nearest neighbors were addressed [3| , finding 
spin-Peierls states, as were suggested from topological events in SU(N) extensions of quantum 
antiferromagnets [j, [^. Dimerized states also result from frustrating longer-range interactions 
as in the case of the Shastry- Sutherland model [Q], where nearest neighbor (J) and next-nearest 
neighbor (J') interactions on alternating plaquettes lead to a VBS for large enough values of J' /J 
{J' / J > 1 in the isotropic case for spin S — 1/2 and in the absence of a magnetic field). In all the 
cases above, a quantum phase transition (QPT) between a Neel state or some other intermediate 
state, and a VBS sets in by varying the ratio of the competing couplings. In particular, in the case 
of the Shastry-Sutherland model, several intermediate phases between the Neel state and the VBS 
were proposed (7l-[Toj. However, the nature of the intermediate phase is still not clear pH - fl^. 

Two-dimensional dimerized quantum antiferromagnets were also discussed recently p^ . [l6j in 
connection to the proposal of deconfined quantum criticality [l7j . In this theory, the critical point 
is characterized by the deconfinement of new degrees of freedom resulting from fractionalization 
of the order parameter, as opposed to the Landau-Ginzburg- Wilson (LGW) paradigm of phase 
transitions, where the critical point is characterized solely by the critical behavior of the order 
parameter and its correlators. Specifically, the theory by Senthil et al. (l7| allows for a continuous 
phase transition between a Neel state, that breaks rotation symmetry in spin space, and a VBS, that 
spontaneously breaks the lattice symmetry, while on the grounds of an LGW-theory, a first order 
phase transition is generally expected. The quantum Monte Carlo simulations in Ref. [l^ show for 
different dimerization patterns of an S' = 1/2 antiferromagnetic Heisenberg model on the square 
lattice that, while most of the dimerization patterns lead to critical exponents consistent with the 
universality class of the three-dimensional 0(3) Heisenberg model, a staggered pattern leads to 
deviations from it. Hence, the authors raised the possibility of deconfined quantum criticality in 
this case. 

A QPT from a Neel to a paramagnetic phase can be approached from the magnetically ordered 
side by introducing spin coherent states [18], that lead in general to an 0(3) non- linear cr-model 
[T^ [20| , as in the treatment of the antiferromagnetic quantum Heisenberg model in the context 
of high temperature superconductivity 21 [, or, alternatively, employing a CP"'^ representation [l5|. 
However, the discussion above, shows that it would be certainly interesting to address the QPT 
from the VBS side in the case of a dimerized system. When the VBS consists of nearest neighbor 
dimers, for spins S = 1/2, the dimers are naturally described by the possible states of two spins, 
i.e. by the manifold SU(2) ® SU(2) ~ S0(4). Some time ago, bond-operators were introduced 
[2a - l24| for that purpose, that offer an intuitive picture for singlets and triplets. The application 
of such a representation with different approximation schemes, mostly based on mean-field th eory 
led to rather good descriptions of spin-liquid states in a variety of situations in frustrated [H, [2^ . 
dimerized [2^j2^ , and bilayer antiferromagnets [275 • The aim of the present work is to formulate a 
path integral dealing directly with the S0(4) manifold on the one hand, and with bond-operators 
on the other hand, enabling a field-theoretic treatment of the transition between dimerized and 
other possible states. 

We consider here two equivalent approaches. The first one deals with the generators of the SO (4) 
algebra. While previous work !28| used coherent states, we perform the path-integral quantization 
using the Faddeev-Jackiw symplectic formalism for quantum field-theories with constraints [29l.l30j. 
The treatment of constrained systems was initiated by Dirac f3ll| and continued by Faddeev and 
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Jackiw [29] . These methods were used before for the Heisenberg and t — J model [32| - |35| showing 
consistent resuhs with other works were coherent states were used [s^ . The formal development is 
shown in Sec. [Ill The second approach, based on a path-integral representation of bond-operators, 
is described in Sec. lIIIl In both cases, special emphasis is given to the determination of the measures 
and constraints appropriate for the respective algebras. We expect in this way to allow for the 
treatment of fluctuations beyond the mean-field approximation, at the same level as in previous 
treatments leading to the 0(3) non- linear cr-model starting from a Neel configuration. Section ITVI 
illustrates how the obtained path integral representations may be used to reach a field-theory for a 
spin-Peierls chain. This example was chosen because the effective field-theory is characterized by 
the presence of a topological term that ensures that in the absence of dimerization, the spin-gap 
closes. In Appendix A we show explicitly, how to transform from one representation into the other. 



II. PATH INTEGRAL FORMULATION FOR SO (4) FIELDS 
A. SO (4) operators for a bond 

We start by considering the four states for a bond joining two 5=1/2 states, 

|o) = -i=(in)-i;t)) , 
|i) = -^(in>-iu)) , 



2) 
3) 



1 

71 



(in>+ \m , 
(in>+ i;t)) . 



We can introduce the Hubbard or X-operators ^3f\ defined as 

with = 0, . . . , 3. Trivially, they obey the following commutation rules: 



and the completeness relation 



(ii.i) 

(11.2) 
(II.3) 
(II.4) 



Since we are considering states of SU(2) (g) SU(2) ~ S0(4), we have to construct the six generators 
of SO (4) out of the 16 operators X^^'^ . These generators can be written as 



(II.5) 



where the Latin indices run over 1,2,3. Using the commutation relations (|II.3p we obtain 



(11.6) 



the commutation relations of the generators of SO (4) [38 
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The generators of SU(2) for each site of the bond can be also constructed in the following way: 

^(2) ^li'T''- S") , (II.7) 

where 1 and 2 denote the two spins making up the bond, such that any spin Hamiltonian can be 
expressed in terms of the generators of SO (4). 

Finally, we consider the Casimir operator in SO (4). The Casimir operator, which commutes with 

the generators, is given by the sum of the squares of the generators [SOj i.e. 7" + <S . From the 

. . - 2 ~ 2 

relations (jll.7p . we have T + S ~ 1. There is another bilinear form of operators that commutes 
with all the generators, namely T ■ S = 0. 



B. Faddeev-Jackiw theory and path integral for SO(4) fields 

In order to develop a path integral representation for quantum system whose operators satisfy 
general Lie algebras, the coherent-state approach Il8| is often used. Here we will instead use the 
formalism introduced by Faddeev and Jackiw (FJ) [29l. Isol . l40l | , that provides a way to obtain a 
classical theory consistent with the algebra of the quantum problem without need of dealing with 
the structure of the corresponding group. In this ap pro ach, no formal distinction is made between 
different forms of constraints like in Dirac's theory [3l|, where primary and secondary, first class 
and second class constraints appear. As shown by Faddeev and Jackiw [21, [13, IS^] constraints can 
be incorporated iteratively (see also this section and Sec. IIIIBI) . until the basic brackets for the 
fields can be determined. Once the classical field-theory is obtained, quantization can proceed via 
a path integral or in a canonical way. 

Our purpose is to construct a classical Lagrangian first order in time derivatives of the SO (4) 
fields. The terms containing time derivatives will give the corresponding Berry phases. Following 
FJ, we call these terms canonical. Since the FJ-formalism is not widely used, we will discuss it in 
some detail, such that the presentation is to a large extent self-contained. Given the Hamiltonian 
H , that depends on fields yA, we start by writing a Lagrangian first order in the velocities i/a, 

C{yA,yA)=Y.yAKA{yA)-H{yA) . (II.8) 

A 

Then, the associated Euler-Lagrange equations of motion are (summation over repeated indices is 
assumed) 

'dKB OKa 



dyA dys 



yB~^ = (IL9) 



dy_ 



A 



If the matrix defined by 



Mas^'-^~'^ (11.10) 



dyA dys 

is nonsingular, it is possible to write eq. pi.9[) as 

f) J-f 

VB = {Mab)-'j— (ILll) 

OVA 

On the other hand, from the Hamiltonian formalism the equations of motion are 

yB^{H,yB}^j^{yA,yB}, (11.12) 
oyA 

where {A,B} denotes the Poisson bracket. Then, comparing pi. Ill) and pi.l2[) . we see that 
(Mab)'^ plays the role of the basic bracket (or generalized bracket) of the Faddeev-Jackiw theory, 
i.e. 

{yA,yB}Fj^ (Mab)-' (11.13) 
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These brackets agree with the Poisson brackets for an unconstrained theory. The generalized 
bracket between two quantities F{yA) and G{yA) is defined by 

From (|II.13|) and (|II.14|) it is straightforward to show that the generahzed brackets verily ah the 
properties ol the usual Poisson brackets. 

In our case, we assume that the first-order Lagrangian can be written in terms ol the S0(4)-fields 
as: 

£{r, S) ^ Alt'' + ytf 5° - , (11.15) 

where and are unknown coefficients which must be determined. Since the S0(4)-fields 
must verily the conditions discussed in Sec. Ill Al we are in the presence of a constrained theory, 
where the potential y'") is 

= i/ + af^i+6f^2 , (11.16) 

with H the proper Hamiltonian, and ^2 Lagrange multipliers, and 

= + 52 - 1 , 

= r -5, (11.17) 

are the invariants discussed in Sec. Ill Al which must be considered here as constraints between the 
S0(4)-fields. 

For the classical Lagrangian (jII.15[) . the set ol classical variables in configuration space is {j/a} = 
{T", 5", ^1, ^2}, the coefficients arc {^^a} = {-^ ^-^d the corresponding equations ol motion 
are 

MAByB = ^ — . (n.l8) 

OIJA 

Alter constructing the matrix Mab , it can be readily seen that it is singular because the coefficients 
Ka are independent of the variables 1^1, ^2- This can be remedied by promoting the constraints 
into the canonical terms [40|. For that purpose, we notice first, that when multiplying pi.lSp by 

eigenvectors v^^ corresponding to the zero modes, we obtain 



V 



dyA 



, (n.l9) 



putting in evidence, that the zero modes ol Mab encode the information of the constraints. For 
consistency the time evolution of the constraints should also obey 

f2, = 0, i,j^l,2. (IL20) 

These conditions can be incorporated into the Lagrangian with Lagrange multipliers Ai^2- In this 
way, discarding total time derivatives, the first iterated Lagrangian results [2^, |4^ 

£(1) (r, S) = Ajf'^ + A^S^ + Al ill + A2 ^2 - V-^^ (IL21) 

where V'-^^ = V'-°^\ O;=o = H. The variables ■^1,^2 have disappeared because the constraints can 
now be imposed on the canonical part. 

The new set ol variables is yA = {T°',S"', Ai, A2} and the new 8x8 matrix Mab can be written 

as 





M 


Fi 






F2 




-FT 






Mab = -M' B^ F2 \ , (IL22) 
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where given by 



Bjb - 



dAT dAT 



dTa dTb 



^ = Eabc (Vr X 



and B given by 



o dAf dAi f„ 



dSa dSb 

are 3x3 antisymmetric matrices. Furthermore, Af is a 3 x 3 matrix given by 



Fi is a 2 X 3 matrix given by 



g^f dAl 
dTa dSb 



2Ti Si 
2T2 S2 
2T3 S3 



(11.23) 



(11.24) 



(11.25) 



(11.26) 



and F2 is as Fi with T and S interchanged. 

The unknown coefficients Aj = Aj{T,S) and Af — ylf (T, «S) must be determined in such a 
way that the matrix Mab results nonsingular and the S0(4) fields verify the following relations 











—i 


f'',S'' 




~i 





(11.27) 

These are the classical version of the commutation relations (|II.6p . From (jII.13p and (jll.27p we 
have 

/ R^^r"}^^ {T^,S'}pj {TaAl}FJ {Ta,X2}Fj\ 
{S°',T''}pj {5°, 5^}^^ {Sa,Xl}pj {Sa,X2}FJ 



(Mab) 



F,J 



{Ai,A2} 



FJ 



\^{Tb,Xi}pj — {Sb,Xi}pj — {Ai,A2} 



(11.28) 



FJ 







where in pi. 281) the elements {Ta, Xi}pj, {Sa, Xi}pj and {Ai, X2}pj are unknown but unnecessary. 
For the determinant of Mab , we obtain : 



detMAB ^ 



where D = Di + D2 with 



Mab {VsxA- 



Vrx^n Mb, 



("5 X T), 



(11.29) 



(11.30) 



and 



Do 



V5 X A^^ ■ 5] [(Vr X A'^^ ■ S 



V5 X A* • r 



(Vr X A'^) • "7"] + 2 ('^"'^^ " ^^^'^^ eabcedefMebMfc . (11.31) 
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Imposing the identities pi. 271) in (jll.281) , and after a long but straightforward algebra it is possible 
to show that the coefficients must satisfy the following set of equations: 

(Vr xA^ + VsxA^)-T+ (Vr x + Vs x A^) ■ S ^ -2 , 
(Vr xA^ + VsxA^)-S + (Vr x A^ + Vs x A^) ■ T = , 

(Vs xA^-VtX A^) ■ (T X 5) + (Alab + Mba) (TaTb - SaSb) 



+ Maa {S' -T')+J2 ^~^aa {S^ - T^) = , 
a 

D = -1 . (11.32) 

A possible solution compatible with the equations above is 

VtxA^ = VsxA^, VsxA^ = VtxA'^, Mab = -Mba, (11.33) 

with 

Vt X = -T , VsxA'^ = -S, Mab = -SabcSc , (11.34) 

as can be easily verified. We notice furthermore that, from the last two equations in (jll.34p we 
have 

M - ^ m 35) 

The solutions displayed above are, however, not the most general ones. In the Appendix we show 
explicit forms of A/ and AP , that obey eqs. (III.32P but not all eqs. (III.34[) . Nevertheless, it is also 
shown in the Appendix, that both forms lead to the same Berry phase in the gradient expansion. 

Having obtained an effective Lagrangian with a set of coefficient s , A^ and constraints (jII.17p 
we can write the path integral for the partition function as usual [4l| 

j VTVS 5{T-S)5 (r^ + 5^-1) e"^ , (11.36) 

with the Euclidian action 

y" dr I ^ (aJ ■ drTj + Af ■ drSj^ + H [r, 5] I . (11.37) 

Here we would like to remark a few features. In the measure of the path integral pi.36p only 
the two constraints enforced by (5-functions appear. No other functional between fields is present 
because the determinant of the matrix Mab is a constant (D — —1). The ffi'st term in the action 
pi.37p . i.e. the Berry phase, depends on the coefficients A'^ and A^ , that have to fulfill (III.32[) . 
Hence, they are defined up to a gauge transformation, 

^r^^r+H, ^f^^f+^, (11.38) 

with A = A(^^, A^) a scalar function, in a similar way as in the SU(2) case. In this way, we 
have obtained a path integral for the 5*0(4) fields that is fully consistent with the algebra of the 
quantum problem. 

III. BOND-FIELDS FORMULATION OF BOND VARIABLES 

A. Bond-fields for spin-states on a bond 

It is usual for constrained systems, where it is necessary to deal with operators that do not sa tisfy 
canonical commutation rules, to introduce slave particles by decoupling the original operators [42|. 
In the frame of the Hubbard operators A^'' defined in (111.21) . the following decoupling can be used, 

X^' = tit, , = sh, , (HI.l) 
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with a, & = 1, 2, 3, and where {t^, tjj, s, s^} satisfy bosonic commutation rules 

Sat, [S,S^]=1, [s^^a]=0. (III.2) 



ta4 



Using (jlll.ip it is easy to show that the commutation rules (|IL3[) are satisfied. In addition, the 
completeness condition (|IL4p can be written as 



s 



+ 4ia = i- (ni.3) 



From eqs. pi.Sp . (111.71) and (IIII.I|) we obtain 

^("2) = I (-S^ta - - ieatctlt,) , (III.4) 



that exactly correspond to the relation between spin- and bond-operators introduced by Sachdev 
and Bhatt UM- 

As the path integral (jll.36p shows, only four real variables are independent. Therefore, on passing 
from bond-operators to bond-fields, four constraints are needed. One of them is the completeness 
condition 

^1 =S*S + Oa-l = 0. (III.5) 

In order to obtain the remaining conditions, we use the well known CP^ representation ^] for the 
spin-fields S'(i) and S(2), 

S^.^^^za'^z, Sf.^^^Coa'^u;, (III.6) 

where cr" are the Pauli matrices. The CP^ fields, z = {z^, Zj) and ui = (wj', Wj), fulfill the conditions 

zz = 1 , WW = 1 . (III.7) 

On the basis of eqs. (jIII.4p and (|III.6p . the following relations between bond-fields and CP^ variables 
can be obtained: 







s = 








h = 




t2 = 


i 








ts = 





(III.. 



These equations lead to two additional constraints for the bond-fields: 



ip2 = SS~ tj^ = , 

(^3 ^ -3*8* ~ tX = . (III.9) 

As can be seen from (|III.4p . the theory defined in terms of the bond-fields contains a gauge degree 
of freedom such that, the remaining constraint will appear as a gauge- fixing condition, that we 
discuss in the next section. 
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B. Faddeev-Jackiw theory and path integral formulation for bond-operators 



In this section we will develop the FJ-formalism for bond-operators. The procedure will parallel 
that of Sec. Ill B[ but due to a gauge freedom, it will be iterated, until the condition for gauge 
fixing is incorporated in the theory. 

Our starting point is the classical Lagrangian 



(III.IO) 



where V^^^ — H{s*, s, t*,t^) + £,i^pi , i/(s*, s, , t^) is the spin Hamiltonian written in terms of the 
bond-fields, S,i are Lagrange multipiers and the constraints (pi with i 1, 2, 3 were defined in eqs. 
(|III.5|) and pil.91) in the last subsection. As in previous bond-operator treatments [1^, we have 
adopted the usual kinetic Lagrangian for bosons. As we will see below, this is consistent with the 
50(4) algebra. 

From the Lagrangian (IIILlOp . the set of classical variables is {ha} = {qa,Ci}j where = 
{s, s*, ij,, t*}, with the index a = 1,...,8. The corresponding matrix Mab is singular, because 
the coefficients Ka do not contain the variables ^i. In this case the matrix Mab has three zero 
eigenvectors . Multiplying the Euler-Lagrange equations (|II.18p by v^^^\ we obtain as in Sec. 

m\ 

^i' ^r— = -iir- = , (ni.ii) 

OVA ail 

that leads to ipi = 0. As described in Sec. Ill Bl and in Refs. [l^, liO] we incorporate the constraints 
()III.11|) into the kinetic part of the Lagrangian using new Lagrange multipliers \i . The first iterated 
Lagrangian results 



(III.12) 



where V'^^'^ = y^^^l^.^o = H. For the Lagrangian (jIII.12[) the new set of variables is ua — {qa,\} 
and the new matrix Mab is 



Mab = 



where 



fafS - 



J 







(III.13) 






i 








—i 




















iSab 








-iSab 






(III.14) 



Note that f^^^ is the matrix formed by the classical extension of the commutations relations for 
bond fields. 

As the matrix Mab is antisymmetric with an odd number of rows and columns (11 x 11), it is 
singular. After the first iteration, Mab has only one zero-eigenvector u [i^ 



(III.15) 



As before, the new constraint is obtained multiplying the Euler-Lagrange equation (|II.18|) by the 
zero-eigenvector u, i.e. the new constraint is 



= ul 



dyA dqa dqp 



(III.16) 



where we have used the equations piI.15|) and (jII.14[) . Therefore, no new constraint arises from 
equation piI.16p because tpi is already a constraint and its time evolution is also a constraint for 
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consistency. This fact is related to the existence of a gauge degree of freedom associated with the 
completeness condition (pi. 

In order to obtain a gauge fixing condition ip^ as a new constraint, we will repeat the process 
introducing a new Lagrangian C^^"^ 

6^\T,S) {h*S - SS* + iX - tat:) + hV: + A4^4 - V^^^ , (III.17) 

where V'^'^^ = F'-^^l^^^o = H and the set of variables of the new configuration space is {ua} = 
{(7a,Ai,A4}. We will choose a gauge fixing condition (^4 in such a way that the new matrix Mab 
is nonsingular. In other words 

n 2 



det [Mas] = 16 



^^4 ^ _ 9^ 
ds ds* 



dta ^ dt* ''^ 



= 16 



(III.18) 



must be different from zero. Computing the inverse of Mab, we obtain the following FJ brackets 
between bond fields: 



{s,s}fj = {s*,s*}fj = , 
{s, s*}fj = i (1 — ss*) + {s wl — s*wi) 

{S^tajpj = taWl ~ SVa , 



{s,tl}FJ 
{s*,ta}FJ 
W,t:}FJ 
{ta, tb}FJ 
{ta,tb}FJ 
{t*a,tl}FJ 



where 



Wl 



is*ta + isV*-t*aWi) , 
-istl + {s*Va - t^wl) , 

t>i - s*v: , 

- (taVb - hVa) , 
i{Sab-tlti,) + {t^V^ -tlVa) 

- {tivt - tiv:) , 



(III.19) 



(III.20) 



and a = 1,2,3. Note that on taking into account the constraints, the FJ brackets are different 
from the usual bosonic commutation rules. As expected, using (IIII.lQp . for any explicit form of a 
gauge fixing, the S0(4) algebra is fulfilled. 

Finally, the partition function for the bond-fields can be written as [4l| 



Vs* VsVtl Vta {det [Mab])^^^ S[ip4 S[ipi] S[ip2] S[ip3] e 



where 



Jat [{s*s + tli,) + H] 



(III.21) 



(III.22) 



and {det \Mab\Y^'^ = 4 A is equivalent to the Faddeev-Popov (FP) determinant IS.fp in gauge 
theories [4l| . 

Before closing this section, we would like to remark that bond-operators were frequently used 
at the mean field level, or closely related approximations, where the full form of (jIII.21[) is not 
respected [23l - [27j . There, only the effective Lagrangian L and 1^1 are consider while tp2, (ps, (p4 
and Afp = {det [Mab])^^^ are missing. While on a mean-field level the measure and constraints 
are not crucial, their presence is important when considering the effect of fluctuations. We show 
in Sec. IIV B 31 how, dealing with them, leads to the correct long- wavelength effective action. 

We will show in the following with the one-dimensional spin-Peierls systems as an example, 
how the corresponding continuum theory can be recovered starting with a path integral for the 
dimerized state. 
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IV. SPIN-PEIERLS CHAIN 



As an application of the formulations developed above, we consider a spin-Peierls chain. On the 
one hand, this is a well known system. On the other hand, the corresponding field-theory has a 
topological term, closely related to the one present in the field-theory for the antiferromagnetic 
Heisenberg model. In order to obtain it, a proper treatment of the fluctuations is needed. 

The spin-Peierls Hamiltonian in one dimension can be written as 

5P = [1 + (-!)'■ A] (IV.l) 

i 

where A < 1 indicates the degree of dimerization. For A = 1, the system breaks down into a set 
of decoupled dimers, while for A = 0, it reduces to the antiferromagnetic Heisenberg chain. 

We can introduce strong bonds given by -I- 1) with i = 2j, j G Z. Using this notation the 
spin-Peierls Hamiltonian (IIV.1[) reads 

Hsp^JY. [(1 + ^) ^.-,(1) • ^.-,(2) + (1 - A) 5,,(2) • , (IV.2) 

i 

where the dimerization pattern is described by bonds labeled by the index j. 



A. Continuum limit with SO (4) fields 



Close to the point where the system goes over to the state appropriate for the antiferromagnetic 
Heisenberg chain, a large correlation length should be expected, such that the continuum limit is 
appropriate. Using pi.7l) . Hsp can be written in terms of the S0(4)-fields as follows 



J 



SP 



E 



(l + A) ij^-Sf, 



(1-A) (r,-5,)-(rj+i + 5,+i) 



According to pi.36|) . the path-integral for this Hamiltonian is given by 
Z = j VTVS 6{T-S) S (r^ + 5^-1) e-^ , 
with the action in imaginary time 

/ dr I ^ (^]" • drTj + Af ■ drSj) + Hsp [T, S] | 



(IV.3) 



(IV.4) 



(IV.5) 



We approach the continuum limit by performing a gradient expansion around a configuration that 
is appropriate for large A, where the Hamiltonian is dominated by 



i/A = ^$:(i + A)(r^5 



(IV.6) 



From a mean-field point of view, the lowest energy is obtained by maximizing and consequently 
minimizing 7"^ • The classical configuration of lowest energy is given by = 1 , so that we take 

(IV.7) 



Sj — CjUj 



with rtj = 1, such that the condition + S"^ = I leads to 



Since Tj is proportional to the change of Sj, i.e. 



(IV.8) 
(IV.9) 
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the fields entering the action are 



(IV.IO) 



where a is the lattice constant of the new lattice. In addition, € • n = 0. 

After defining the fields for the gradient expansion, we consider the different pieces of the action. 
First we have 



(IV. 11) 



Hence, the modes described by the field £j are massive. Since we are in 1+1 dimensions, we perform 
the gradient expansion up to O (a^). For the rest of the terms coming from the Hamiltonian we 
have. 





Tj+i 








■ ^j+i 


~ a^£j ■ dxTij 


1 




Tj+i 








■ ^j+i 


~ 1 - a^l] + 





Putting all the contributions together, and going over to the continuum, we have 
HsP [T, S]^Jaj dx \e + (1 - A) • d^n + ^ {d^nf 
Next we deal with the Berry phase. Here we have 

sb=-i J^^j2 (-^r ■ ^-"^j + ■ ^-^^ 



(IV.12) 



(IV.13) 



(IV.14) 



Since from the discussion above, the S0(4)-fields have only a smooth spatial dependence, we can 
concentrate on one site. Furthermore, since from (IIV.IOI) | 7" |~ a, we can expand the vector 
potentials accordingly. In lowest order we have from pi. 341) . 



Vr X A 



T 



r=o 



r=o 



(IV.15) 



Hence, we can write 



iT 



(r,5) 



d 



(r = o,5) = — 0^(r,5) , 

where (ff^ and (f)^ are scalar functions. Such terms can be gauged away when necessary. 
The expansion of the fields A^ and A^ in powers of T° leads to 



(IV.16) 



{T,S) = At iT=0,S) 



dAf, 



r=o 



(IV. 17) 



r=o 



where due to the presence of one time derivative in the Berry phase, we need to consider only 
terms up to 0{a). Introducing the expansion into (jIV.14p leads to 



A^ ■ drT + A'' ■ drS = A' (0, S) ■ drT + 



T , 



dAT 



+A^ io,s)-drS + 



r=o 



(IV.18) 



r=o 



13 



From PL34I) we have 



^•^r _ ^-^T I abcq-c 

Qj-b dT"- ' 



(IV.19) 



Qq-b Qq-a 
such that using (jll.SSp . we arrive at 

■ drT + ■ drS ^ dr (a^ -r^+A^ (0, 5) • drS + e''^''drS''T''S'' 

= e^'^^a^^'^r^^" , (IV.20) 

where the first term can be discarded since it is a total time derivative and the second one can be 
gauged away, as discussed after eq. PV.16|) . Inserting the relations pV.10[) . we finally have for the 
Berry phase in the continuum limit, 

SB=-i j drdx^ • (n X drn) . (IV.21) 

Going back to (|IV.5p . we have for the action after the gradient expansion. 



5 = - 



j drdx^ • (n X drn) 



+Ja / drdx + {1 



1 1 2 
-i ■ dxn H — (dxTi) 

2 8 



(IV.22) 



At this point £ can be integrated out. Since the action is quadratic in this field, we can simply 
consider the saddle point for {., that leads to 



1 

2Ja 



Ja 



(1- A)a,n'^ 



(IV.23) 



Inserting this into the action, we finally obtain 
S = [drdx 

.(1-A) 



+1 ■ 



drdxn • x d^n) 



The coupling constant and spin-wave velocity are given by 

4 

g 



(1 _ A2)l/2 



and 



2 L 



.211/2 



(IV.24) 



(IV.25) 



(IV.26) 



The same result can be achieved performing the same treatment as for the Heisenberg model |42| . 

The action (jIV.24|) corresponds to an 0(3) non-linear cr-model with a topological term i9Q, with 
9 = S{1 — A)/2, for 5* = 1/2, and Q given by the integral in the last term of (jIV.241) . acquiring 
values Ann, with n an integer. For the antiferromagnetic Heisenberg chain (A = 0), it is responsible 
for the absence of a gap for half-integer S, while a gap is present for S integer, as first discussed 
by Haldane [3, some time ago. Affleck and Haldane [46l - l48j related the 0(3) non-linear cr- 
model with a topological term to the Wess-Zumino-Witten non-linear cr-model with topological 
coupling constant A: = 1, a conformal field theory that constitutes an attractive fixed point for 
antiferromagnetic Heisenberg chains with half-integer spin. Based on the conformal field-theory, 
the critical behavior for A 7^ was discussed by Affleck et al. [49], showing that the mass gap 
opens as m c>c A2/'^/| In Aj^/^^ Numerical simulations with a cluster algorithm 'so'] confirmed later 
the correctness of the identification of the action pV.24l) with the WZW conformal field theory, 
and showed that for A 7^ 0, a gap opens. Hence, our path integral formulation allowed us to 
start from a dimerized phase and take into account the relevant long-wavelength fluctuations that 
describe the critical behavior of the system in the limit A — >■ 0. 
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Continuum limit with bond-fields 



As a first step, we have to identify as in tfie previous discussion, slow and fast components of 
the bond-fields, in order to allow for a gradient expansion. For that purpose, we use a staggered 
CP^ representation, that was previously introduced for the t-J model [3a|. We discuss it shortly 
here, in order to allow for a self-contained presentation. 



1. Staggered CP^ representation 



We consider now the fact that the spin-fields S'(i) and S'(2) on a bond are staggered with respect 
to each other, and introduce a CP^ representation, such that 



Qa 



^(2) a a V (2)* 



where z*-'-* z*-*-* = 1, i = 1, 2. On the other hand, making the following replacements 



.(2) 



Z2 
UJ2 



Z2 



in the expressions (jIII.8[) . the bond- fields can be written as follows: 



2 

,^^(2)^a^(l) 



From the relations above it is easy to show that the completeness relation 



s*s + t*„t„ 



,(1)* (1)^(2)* (2) 



and that the constraints 



also hold. 

The Berry phase is given by 

This is the form of the Berry phase in the staggered CP^ representation 



(IV.27) 



(IV.28) 



(IV.29) 

(IV.30) 
(IV.31) 

(IV.32) 



2. Slow and fast components of the bond-fields 
For each bond we can define new fields 



aC 



Due to the constraints satisfied by the fields z'-'-', we have 

(I -I- aC) {i + aQ^l , 
{z - aC) (z - aC) = 1 , 



(IV.33) 



(IV.34) 
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leading to constraints for the fields z and C, 



Iz + a^CC = 1 • 

We then introduce new fields 

z = z^/f - a^CC , 
such that the constraint (|IV.36|) goes over into 

zz = 1 , 

and the constraint pV.35|) translates into 

zC + C^ = . 

The fields z and C defined above correspond to smooth configurations. 

Next, we can express the bond-fields in terms of the fields z and Q as follows: 

s = [1 + a (Cz - zC) - 2c?lC\ , 



(IV.35) 
(IV.36) 



(IV.37) 
(IV.38) 
(IV.39) 



/2 

= il-^ [za^z (1 - a^CC) + a (Cct°z - zct°C) - a^C^^C] 



(IV.40) 



where we keep contributions up to O (a^). Guided by the expression of ta, we define the following 
vectors 

O = zcrz , 
L = zcrC , 
= Co-z , 
m = C,cr( . 



For later convenience, we define the following two vectors in addition to f2 and m 



1 



Ln^-iL + L') =-{zctC + Cct 



1 



L L 



Further relations can be obtained by using that 

ft- L^zC, n- l'' ^(z , 

such that 



(z - zC ^ n ■ \^L' ~Lj = ~2iQ. ■ Li . 
Using the constraint pV.39[) it can be also shown that 

With the relations above, we can rewrite (jIV.401) as follows 

s = [l - 2ian ■ Li - 2c? {Lif \ , 

ta = [^a - 2iaLia - (i/fia + "1q)] 



(IV.41) 

(IV.42) 
(IV.43) 

(IV.44) 

(IV.45) 

(IV.46) 



(IV.47) 
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Due to the constraints (jIV.38|) and (jIV.39|) . the vector-fields fulfill the following conditions: 







= 1, 


n 


Lb. 


= 0, 


Li 


Lr 


= 0, 


m 


Lr 


= . 



(IV.48) 

There we can see that the vector-fields J7, Li, and m are coplanar, so that they can be described 
as linear combinations of two orthogonal fields in that plane. We choose them to be O and O x Lr. 
Using the form of r2 x Lr expressed in terms of the fields z and ( (eq. pV.41l) and eq. (|IV.42p ) 
we can write 



Li ^{n- Li)n- (ft X Li 



and 



m = 



2 {ft ■ Li f -L] rj - 2 (O • i/) (O X Lr) 
Inserting the expressions for Li and m into (IIV.47|) . we have 



(IV.49) 
(IV.50) 



n/2 



[rt - 2m [ipQ,- [rtx Lr)] - 2a^if [(^$7 - (fJ x Lr)]] 



(IV.51) 



where we introduced Lp = Vl ■ Li. Here we see that using the fast and slow fields of the staggered 
CP^ representation, we obtain an expansion of the bond-fields in terms of two vector- and one 
scalar-field. The constraints PV.30|) and (jIV.31|) on the bond-fields are fulfilled by imposing the 
conditions (|IV.48|) . such that there is no condition on the scalar field at this stage. 
Recalling that the bond-fields are related to the S0(4) fields as follows: 



Sa — s*t^ -\- t*^s , 

we can relate the SO (4) vector- fields with the fields entering (|IV.5ip . arriving at 



(IV.52) 



S 

r 



n - 2a^ 
-2aLR 



L%^ft + (p{n X Lr) 



(IV.53) 



One can easily see that the fields above fulfill the constraints obeyed by the fields S and T up to 



3. Constraints and measure 



The expansion (|IV.5ip of the bond-fields in fast and slow components involves only seven com- 
ponents instead of eight, as originally introduced. Taking into account the first two constraints in 
(|IV.48p , we are left with five fields instead of the four required to describe the degrees of freedom 
on a dimer formed by two 5' = 1/2 spins. Therefore, before passing to a path integral in ip, fl, and 
Lr, a gauge fixing condition is still necessary. Here we choose the gauge fixing ^3 -I- ^3 =0. Once 
we proposed an expansion of the bond- fields in terms of fast and slow variables, we consider first 
the constraints on those variables together with the change of measure due to the transformation 
to these new variables. 

The measure and constraints corresponding to the general gauge-fixing were obtained in Sec. 
HHBl fcq. IIII.2ip . For the particular gauge — + ^3 the measure can be written as follows: 



VM - Vs* Vs Vtl Vt, (t* - tg) 5 {tl + 1,) 

X<5 (S*S + tX - 1) S {S*S* - tUl) 6 {SS - tata) 



(IV.54) 
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Since we have in total seven variables given by (/s, $7, and L^, we can eliminate one of the variables 
by integrating over ^3 and imposing the constraint due to gauge fixing. Then, we have 

VM Vs Vtl Vt^ <3 5{s*s + tlt^ + tm -tl-l) 

x5{s*s* ~ tlt^ - t^t^ - tfj 5 (ss - tata) • (IV.55) 

The arguments of the ^-functions with the new variables look as follows 



- 1 



Rl) 



5{i 



W - 1 



2ia {{l-n')^ + [^n-s - {^iLr2 - f^s^m)] } 



+2a2 

5 [ss - tata) 



R) 



(n^ - 1) (Ll 



- 1 



2iaip - 1) + 2a^ (J7 • Laf - {l\ + 2ip^) {fl^ - l] 



We can transform these constraints onto constraints on the new fields: 

6 {s*S + tlt^ + t*2t2 - ^3 - 1) 6 {s*S* - tlt^ - t^t^ - tg) 5 {ss - tata) 



'J 5 



ifQs - {QiLr2 - ^2Lb.i) S (fl-LR) 
where the inverse of the Jacobian for the transformation is up to O (a^) 



T-i , ■ 2 

^ 2n 



if {2lp^z + fliLB2 — ^2Lm) 



2nl 



(IV.56) 



(IV.57) 



(IV.58) 



'3 ^"3 

On the other hand, the Jacobian for the transformation to the new fields is up to O (a^) 
J = I6V2 (O • Lr) |i73 - 2ia [cpQ^ + (f^ii^ja - ^^2^^)] 

-2aV [^^3 + {^iLr2 - ^2Lri)] } . (IV.59) 

The first factor is cancelled in going from 6 Lr)^ io 5 {VL- Lr). Finally, the measure in 
PV.SSP is now 

{D,iLr2 — D,2Lri) 



VM cx Vip Vn VLr 5[W -\)5 



6in-LR). 



(IV.60) 



4. Field-theory for the spin-Peierls model 

Changing notation Lr L, for brevity, we just express the S0(4)-fields in terms of ip, ft, and 
L using (IIV.52I) . where keeping terms up to O (a^) leads to PV.53[) . Going back to (lIV.Sp . we can 
express the Hamiltonian with the new fields, leaving aside constant terms. 



Hsp = I dx'^ 



mL^ + (i-A){ld.n.d..n-4d.n.L 



(IV.6I) 
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Next we consider the Berry- phase as given by (jIV.321) , where we need to consider only contribution 
up to 0{a), and discard total time derivatives, 

Sb= J dTY^{s*s + t*Ja) = -2i J dTdxL-{nx drn) . (IV.62) 
Finally, the total action is given by 



S ^ J dTda;| - 2iL ■ (O x 9^0) 



J a 



16i^ + (1 - A) - d^n ■ d^n -Ad^n ■ L 



with the partition function 

J VipVnVL5{VL^ - I) 5{n-L) 5 
Therefore, we can trivially integrate out ip and obtain 



(riiLa-fiaii) 



Z = JvnVL5{n^~l) 6{n-L) e"'^ , 



(IV.63) 



(IV.64) 



(IV.65) 



with the action given by pV.GSp . As previously done, we can integrate out L, or what is equivalent, 
we solve for the saddle point for L as done in subsec. IIV A[ arriving at 



4Jia 



(IV.66) 



leading to the effective action 

S — drdx 



(1 - A) drdxfl- (drfl x dM) 



(IV.67) 



that, as expected, coincides with the action pV.24l) . 



V. SUMMARY 

We have presented path integral formulations for dimerized quantum antiferromagnets both for 
the S0(4)-algebra obeyed by the operators on a bond, as well as for bond-operators [l^. We 
used the formalism introduced by Faddeev and Jackiw for the quantization of constrained systems 
[29I I30I , I40I I43I , where the presence of constraints is revealed by singular modes of the symplectic 
matrix resulting from the canonical form in the Lagrangian. As opposed to the Dirac treatment 
of constrained systems [31j, no distinction among primary and secondary, first class or second 
class constraints is needed, but an iterative incorporation of constraints determined by the singular 
modes, into the canonical form. For the SO(4)-formulation, the FJ-formalism ensures that the basic 
brackets of the fields corresponds to those of the generators of the algebra. On the other hand, bond- 
operators posses a gauge degree of freedom, such that depending on the form of the gauge-fixing, 
different forms for the basic brackets may result, a fact that is also known from the quantization 
of slave-particle formulations of the t-J model [51|. As an application for both formulations, we 
considered the spin-Peierls model in one dimension, where we obtained the corresponding field- 
theories. This model, albeit simple and well known, was chosen since the gapless phase is realized 
due to the presence of a topological term. The present treatment shows how, by taking into 
account properly fluctuations in a gradient expansion, this non-trivial information is encoded in 
the fields appropriate for dimerized states. Since the path integral formulation presented here is 
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not restricted by dimensionality, we expect this treatment to be useful in dealing with dimerized 
quantum antiferromagnets in higher dimension, as those mentioned in the Introduction. 
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Appendix A: Mapping between the path integral formulations for SO(4) and bond fields 

In the present appendix, we obtain the path integral representation for the SO (4) generators 
(eqs. (III.36|) and (jll.37|) ) from the path integral representation for the bond fields . The starting 
point is eq. (|III.2ip 



Vs* Vs Vtl Vta {det [Mab])^^^ S[^i] S[^i] S[ip2] S^s] e 



(Al) 



where the action is given by eq. (jIII.22[) . with the constraints (jlll.Sp . (IIII.9I) . and the gauge fixing 
condition (p4 = + 1^, as in Sec. (jlVB 3p . Then, from the equation (|III.18p we have 



Next, we introduced in (jA.l[) the S0(4) fields given by (|IV.52p through the identity 

1 = J vsi VS2 VS3 vTi VT2 VT3 V5[% + i{ti t2 - 1; h)] , 

where P is a product of the remaining (5-functions that relate the S0(4) with bond-fields: 

V = d[Si ~ is*h + tl s)] 6[S2 - is*t2 + t; s)] 6[S3 - [sHs + s)] 
X S[Ti + i{t*2 h - ^3 ^2)] 5[T2 + i{tl h ~ tl ts)] . 

It is possible to show that 



(A.2) 



(A3) 



(A.4) 



r S[ss-t^ t„] S[s*s* t*J 



J 



x6 



53 T2 Si ~ Tl S2 



, , Sz T2S1 —T1S2 



t; + 



x5 



t3 Si - i s* T2 
S3 

t3S2-isTi 

S3 



t3Si~ isT2 

S3 



t% + 



t3 ^2 - i S* Tl 
53 



S3{Sl+Ti + T?) 



\ ^{Si Tl + 52 72? + 5| (5^ + + Ti) 



(A.5) 



where J = — 4^3 (^3 — is the Jacobian of the change of variables. Using (IA.5I) . and after 
integrating over the bond-fields, we arrive at 

Z = JvTVSSiT-S)S{T^ + S^-l) e"^ , (A.6) 

where the Berry phase is given by 
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resulting in the following contributions after simplifying the expressions using the constraints of 
the SO (4) fields, 





— 


T2T3 


(Si + Ti 


+m 


Al = 


~SiS3 4 




isi + V 


+ V) 








Al = 


{T2S1 - 


T1S2) 


{Si + Ti 





(A.8) 

where the curl of the fields above only partially fulfill eqs. (jIL33l) and (jll.34[) . However, they 
comply with eqs. pi.32[) . Using the expressions (IA.8I) . and the expansions of the fields S and 7" 
from eq. (|IV.10p . we arrive at the following expression for the Berry phase: 

■ drT + ■ drS = -i— (iidrn2 + lidrn2) + O (a^) . (A.9) 

It is straightforward to show using the fact that — 1 and n ■ £ ^ 0, that the expression above 
is equal to the integrand in eq. pV.21[) . 
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